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^ ! Abstract 

■ Formulas for the contribution of the conduction electrons to the polarization and magneti- 

zation are derived for disordered systems and within a one-particle framework. These results 
generalize known formulas for Bloch electrons and the presented proofs considerably simplify 
Q-f and strengthen prior justifications. The new formulas show that orbital polarization and magne- 

tization are of geometric nature. This leads to quantization for a periodically driven Piezo effect 
as well as the derivative of the magnetization w.r.t. the chemical potential. It is also shown 
how the latter is connected to boundary currents in Chern insulators. The main technical tools 
in the proofs are an adaption of Nenciu's super-adiabatic theory to C*-dynamical systems and 
Bellissard's Ito derivatives w.r.t. the magnetic field. 
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1 Introduction 



Polarization and magnetization are two important properties of solid state systems. In classical 
electromagnetism they are often described phenomenologically, but within the framework of quantum 
mechanics they ought to be calculated from first principle. It turns out that both are composed of two 
main contributions, a relatively simple one resulting from the ions and spins respectively and a more 
delicate one stemming from the conduction electrons. This latter contribution is called orbital and 
this work is about how to calculate it, at least in principle, for disordered solids within a one-particle 
model for which the interactions between electrons (or quasiparticles) are neglected. 

Let us first focus on the polarization and the Piezo effect as the most prominent example. Defor- 
mations of a crystal can lead to a polarization of the sample, i.e. to the accumulation of charges on 
the boundary of the sample. The first contribution to this effect is the polarization due to relative 
displacements of the ionic cores in a unit cell of the crystal and it is straightforward to compute. The 
second contribution comes from the conduction electrons. This contribution was only understood 
in the last 20 years. The key observations are, first of all, that the polarization of the sample is a 
bulk effect resulting from a non-vanishing current induced by a change of the Hamiltonian in time 
and, second of all, that only changes in the polarization are well denned. Based on a linear response 
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argument, King-Smith and Vanderbilt derived a formula for the polarization in terms of instanta- 
neous Bloch eigenfunctions of the underlying periodic system |KS V] . A review of the more recent 
physics literature is given by Resta |Res] . The King-Smith and Vanderbilt formula was rigorously 
derived and generalized to continuous and periodic systems in |PST] . More precisely, it was shown 
that a certain generalization is correct in the adiabatic limit of slow deformations without power law 
corrections. 

In this work, the focus is on independent electrons in a disordered medium and in a tight-binding 
approximation. The class of space-homogeneous Hamiltonians is described in detail in Section 12.11 
Motivated by the above picture, let us suppose that there is a continuous path t G [0,T] h» H{t) of 
time-dependent Hamiltonians. Then states evolve according to the Liouville equation, see Section [2~T1 
The change AP in polarization during a change of the system between time t — and t = T is then 
defined by integrating the induced current density carried in the time-evolved state p(t) of the system: 

ap = f ' dtr B (p(t)j(t)). (i) 

Jo 

Here J(t) = i[H(t),X] is the current operator given by the commutator of H(t) and the position 
operator and 7b is the trace per unit volume at fixed magnetic field B. All these objects obtain a 
clear mathematical meaning in Section [2j Note that AP is still a density and needs to be integrated 
over a suitable surface to compute the accumulated surface charge. Let us also point out that, even 
if the initial condition p(0) is given by a function of the Hamiltonian, it is not true that also p(t) is a 
function of the Hamiltonian (in which case the integrand in (Op) would vanish by Proposition [3] below). 
It is now relatively easy (see Proposition Hj) to derive that, if p(0) = Po(0) is the Fermi projection on 
all states below the Fermi level p(0) lying in a gap of H(0), then p(t) is also a projection and 

AP = [ T dtT B (p(t)[d t p(t),[X,p(t)]]) . (2) 

JO 

This formula is not very useful because it requires the knowledge of the time evolution of p(t) which 
is, in general, not a function of H(t). However, in the adiabatic limit of time slow changes, the 
situation is simpler as proved in our first main result. 

Theorem (detailed version in Theorem [1] below) Let t (->■ H(t) be smooth with a gap that remains 
open in t and let 4zH{t) be compactly supported in (0, T). The spectral projection below this gap is 
denoted by Pq(£). For the time- evolution governed by ed t p(t) = i[p(t), H(t)] with initial condition 
p(0) = Po(0) one has for any N G N 

AP = f dtr B (P (t)[d t P {t),[X,P Q {t)))) + 0{e N ). (3) 
Jo 

The leading order expression in ([3]) is the above mentioned King-Smith and Vanderbilt formula. 
As already pointed out in the abstract, apart from (j2J) the key ingredients to the proof is an adaption 
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of Nenciu's superadiabatic expansions to the C*-algebra of homogeneous observables. In the Corol- 
laries [2] and [3] it is shown that the r.h.s. of ([3]) is of topological nature. It does not change under 
diffeotopic deformations of the path t h-> Po(t) and, for a periodically driven system of period T, it 
is qunatized, as noted by Thouless [Thoj in a more restricted context. 

Next let us come to the magnetization of a solid. First of all, there is a contribution due to 
alignment of nuclear spins. This phenomena can be studied by a Heisenberg model. It may lead to 
a local magnetic field, which in turn can generate a second so-called orbital contribution resulting 
from the motion of the conduction electrons by the Biot-Savart law. If the electron Hamiltonian is 
time-reversal invariant, then the orbital magnetization vanishes (see Section HJ). If, however, there is 
a spin magnetization, the time-reversal invariance of the electron system is broken possibly resulting 
in orbital magnetization. Moreover, it may then be of comparable size, even though a secondary 
effect. Nevertheless, only quite recently there has been progress on how to calculate the orbital 
magnetization in the physics community. One reason for this may be that the magnetization cannot 
be calculated by linear response theory. The first (independent) works seem to be |CTVR[ IXYFN] 
and they exhibited the orbital magnetization as a geometric quantity which can be expressed in terms 
of the Berry curvature and the Rammal- Wilkinson tensor of the Bloch bands. Since then there were 
numerous contributions which are well-documented in three recent reviews [Rest IXCNj IThoj . 

Now let us define the orbital magnetization of a system of independent electrons described again 
by a (covariant family of) one-particle Hamiltonian H. One starts with the pressure (free energy per 
unit volume) at inverse temperature /3 and chemical potential \x: 

PGM = ^T B (ln(l + e-^))) . (4) 

Then the components of the magnetization M(/3,/i) = (M 1 (f3, //), . . . , M d ((3, //)) are its derivatives 
w.r.t. the magnetic field: 

Mj(J3,ii) = d Bj p{^p). (5) 

Of course, one has to verify that the pressure is indeed differentiable. This is obvious within the 
formalism presented below, but can seem rather puzzling if one has in mind that a constant magnetic 
field B asks for an unbounded vector potential for which there is little hope to control perturbation 
theory. The same technical problem appears when trying to calculate the magnetic susceptibility, as 
explained in [BCS] . One way out is to use periodic magnetic fields and then study the long wave 
length limit [MLl ISVXN] . Here we follow a different path and view the magnetic field as parameter 
encoding the algebraic structure, more precisely, the type of non-commutativity of the magnetic 
translations. Therefore, the Hamiltonian is a single element in a field of algebras indexed by the 
magnetic field and within this field derivatives w.r.t. the magnetic field are readily defined. This 
idea was used by Bellissard [Bell] to prove Streda's formula [Str] as well as Lipshitz continuity of 
the boundaries of the Hofstadter butterfly [Bel3j. A new key technical supplement to this formalism 
is a generalized DuHamel formula for the derivation w.r.t. the magnetic field (see Section [23]) . 

The first main result on the orbital magnetization (Theorem [2]) is that Mj(j3, ft) can be expressed 
in a simple manner and without any approximations in terms of the current-current correlation 
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measure (which also appears in the Kubo formula for the electrical conductivity \KP\ ISBB] ) . From 
this we deduce the following new formula for the orbital magnetization. 

Theorem (detailed version in Theorem [3] below) Let the Fermi level // lie in a spectral region of 
dynamical (Anderson) localization and let P^ denote the Fermi projection. If the indices j are under- 
stood cyclically in space dimension d = 3, then the zero-temperature orbital magnetization is given 
by 

M^oo.ai) = ~ M\pi-H\ [[X j+1 ,P^[X j+2 ,P„]]) . (6) 

For Bloch electrons and fi in a gap, it is shown in Section [5] that this formula reduces to that 
of |CTVR] . However, our compact formula is new and we believe the stated hypothesis to be 
necessary. Let us also stress that the presented derivation of (EJ) is considerably more translucent 
than prior ones (apart from being rigorous and more general). The expression on the r.h.s. of fl6]) is 
somewhat similar to the formula for the Chern number in the quantum Hall effect [BESj . Actually, 
the derivative c^Mj (oo, ji) is up to a numerical factor precisely equal to the Chern number as long as 
the Fermi level is in a region of dynamical localization (Theorem H]). In particular, this derivative is 
topologically quantized to an integer and thus constant in any given interval of dynamical localization. 
This shows that the magnetization is indeed of geometric nature. 

In |CTVR] it is already pointed out that boundary currents contribute to the orbital magnetiza- 
tion. Indeed, let us consider a two-dimensional sample in the form of a disc of radius r which has 
a boundary current I flowing around it. Then the associated magnetic moment is of order Ir 2 , so 
that the magnetization (magnetic moment per unit volume) is proportional to I. In Theorem [8] it 
is stated that such boundary currents calculated from first principle in a system on a half-space are 
topologically quantized with the Chern number given by the bulk system. Therefore the strength 
of boundary currents is dictated by (the /i-derivative of) the magnetization itself or inversely the 
magnetization can be read off the boundary currents. This is in the spirit of topological insulators 
and similar to the situation in the quantum Hall effect [KRS] . 

Acknoledgements: H. S.-B. thanks Baptiste Savoie for bringing the work of Resta |Resj to his 
attention. S. T. thanks Gianluca Panati, Giuseppe DeNittis and Max Lein for numerous discussions 
on polarization and magnetization in the past. This work was supported by the DFG. 

2 Formalism 

In order to make this work essentially self-contained, we review the C*-algebraic description of inde- 
pendent particles in homogeneous media [Bell tlBel2] . Even though a bit lengthy, this is indispensable 
in order to fix notations. As already explained above, it is also a convenient framework for the study 
of variations in the magnetic field |Bel2, Bel3j. 
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2.1 The Hamiltonian 



The aim is to model electrons or holes (quasiparticles) in a solid which can be either periodic (so that 
Bloch-Floquet theory applies), quasiperiodic or random (as in the Anderson model). Let us suppose 
that the particles have a spin s and are submitted to a periodic, quasiperiodic or random potential 
as well as a uniform magnetic field. There may also be spin-orbit interactions. The centers of the 
elementary cells of the underlying crystal will be identified with Z d . We restrict ourselves to the 
tight-binding representation describing the low energy states. Hence the one-particle Hilbert space 
describing the quantum states is % = £ 2 (Z d ) £g> C L where L = (2s + l)r and r is the dimension of 
the other internal degrees of freedom (e.g., r = 2 for the honeycomb lattice in d = 2). The main 
focus will be on the cases d = 2 and d = 3. The constant external magnetic field is given by an 
antisymmetric real matrix B = (-Bjj)ij=i...d, namely B^j = —Bj^. The magnetic translations 

(£/aV0n = e^ Bn Vn-a , (7) 

define a projective unitary representation of the translation group Z* d on %. Here the notation 
a ■ Bn = a iBi,i n j is used. For dimension d = 2 and d = 3, the magnetic field is parametrized as 
follows 
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In dimension d = 2, the magnetic flux trough a unit cell is B 3 , and, for d = 3, it is B ly B 2 and B 3 in 
the directions 1, 2 and 3 respectively. The exponential factors in §J\ are explicitly given by 

e fa-Bn _ e %B3(aim-a,2Tli) g|o-Bn _ ^\B\(ainz-azn2) ^\Bi(azn\-a\nz) ^\Bz(a\ri-i-a-i,ri\) ^ 

As the magnetic field B only enters all constructions below through the magnetic translations, it is 
parametrized by values in the torus 5 = [— ir, ix) 2 d ( d_1 ). 

The class of Hamiltonians of interest are strongly continuous families (H^uen of self-adjoint 
bounded operators on % indexed by a compact space Q of disorder or crystaline configurations. This 
space is supposed to be compact and furnished with a continuous action T of the translation group 
Z d . The family of Hamiltonians is then supposed to satisfy the covariance relation 

UaH^U- 1 = H T a u , a e Z d . 

In a typical situation, H w = H Q + is the sum of a translation invariant kinetic operator H and a 
potential energy of the form 

H w = Ho + Vu, H = J2 T mU m , V w = Y,V U)n \n)(n\. (10) 

Here the sum over m has only finitely many terms and the coefficient matrices T m G Mat(L, C) 
satisfy (T m )* = T_ m and (U m ip) n = e~i m ' Bn "ijj n _ m . Furthermore, the matrices V^^ G Mat(L,C) are 
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self-adjoint and continuous in u, and we used the bra-ket notations to denote by \n) the partial 
isometry in H onto L-dimensional state space supported by n G r L d . Then, if the V^ iH are drawn from 
a compact set K C Mat(L, C), the set Q is simply the Tychonov space K xZd furnished with the shift 
action. It is also possible to construct Q as the hull of a single given Hamiltonian H on % |Bell] , 
namely as the strong closure of {UaHU^ 1 \ a G Then the Hamiltonian H is called homogeneous 
whenever the hull is compact. 

2.2 The C*-algebra A of observables 

Now let us construct an abstract observable C*-algebra A the representations of which contain all 
the operators H^, but also functions of H u as well as commutators of with the position operators 
and many more interesting operators which are space-homogeneous. All magnetic fields B varying 
in the torus 5 appear as a parameter in this algebra which allows to study variations in B. The 
C*-algebra A = (*4.B)BeH is continuously fibered by the magnetic field, with each fiber Ab simply 
being a (reduced) twisted crossed-product of the action of Z d on Q. One advantage of this approach 
is that physically interesting objects can compactly be written out with derivatives and integrals 
that are automatically well-defined. In the case of periodic solid, all formulas reduce to standard 
textbook expressions of Bloch theory. 

The construction of Ab goes as follows. One endows the topological vector space C C (Q x 
Z d ,Mat(L, C)) of continuous functions with compact support on Q x Z d and values in Mat(L,C) 
with a *-algebra structure: 

AB(u,n) = A{u,l)B{T~ l uj,n-l)e> m , A*{u,n) = A{T~ n u,-n)* , (11) 

where the last * denotes the adjoint matrix. For u> G £1, a representation of this *-algebra on % is 
given by 

( W^) n = A (T' n oo, I - n) e^ Bn , ^ G f(Z d ) ® C L . (12) 

zez d 

The representations 7Tb,lj are related by the covariance relation 

U a 7TB,uj{A)U a = TTb,T°-lo (A) , a G Z d , (13) 

and are strongly continuous in B and u. Now ||A|| = sup^ e£7 ||vr BLlJ ( J 4)|| defines a C*-norm on 
C c (tt x Z d , Mat(L, C)) and the observable C*-algebra A B = C(Q) x 7L d ® Mat(L, C) is defined as the 
completion of C C (Q x Z d , Mat(L, C)) under this norm. 

Now let A be the continuous field (^4b)b£H of C*-algebras in the sense of Dixmier |Dixj . The 
compactly supported continuous matrix-valued functions on H x Q x Z d are dense in the C*-algebra 
A and their product and adjoint are then given fiberwise as in fTTTT) : 

AB(B,w,n) = ^(B,w,/) J B(B,T-^,n-/)e5"- B ' , A*(B,u,n) = A(B, T~ n u, —n)* . (14) 
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There exists a *-homomorphism pb '■ A — > Ab densely denned by evaluation at B and we use the 
notation tib,ui for 7rB,aj ° Pb • 

Finally, let us exhibit explicitly an element H G A the representations 7Tb,lj{H) of which are equal 
to the operator H u given in ( TTOl) . It is given by the following continuous function on H x Q x Z d : 

H(B,co,n) = T n + • (15) 

2.3 The derivations on A 

On the C*-algebra A exists a <i-parameter group k G T d i— >■ p& of ^-automorphisms defined by 

(p fe A)(B,u;,n) = e lfc " n A(B, w, n) , 

where k ■ n = Yl d j=i^j n j^ namely pk is linear and satisfies pk(AB) = pk(A)pk(B), Pk{A*) = Pk{A)* 
and the group property pkPw = Pk+w- As every ^-automorphism, pk also conserves the C*-norm. 
Its generators V = (Vi, . . . , V^) are unbounded, but closed ^-derivations with domain C^(A) which 
satisfy the Leibniz rule 

V(AB) = (VA)B + A{VB) , A, Be C*(A) . 

The index s indicates that only existence of space- wise derivatives is requested in C^(A) (derivatives 
w.r.t. the magnetic field are dealt with further down). Explicitly the generators are given by 

VjA(B,u,n) = mj A(B,u,n) , A e C](A) . (16) 

Note that these derivations are also defined on each fiber Ab- Furthermore, the sub-algebra C^(A) C 
A is endowed with the graph norm ||A||ci = ||A|| + ^ - ||VjA|| which makes it a Banach *-algebra 
as can readily deduced from the Leibniz rule. If the position operator X = (X±, . . . ,Xd) on % is 
defined as usual by 

(Xj1p) n = Ujlpn , V = (4>n)n&<i e H , 

then 

7TB j A) = i [ttb ,u (A) , Xj] . (17) 

The following result is proved in [Sak, Theorem 3.3.7], but for a slightly more regular function, it 
also follows from the proof of Proposition [2] below. 

Proposition 1 Let H = H* be an element of C^(A). Then for any f G C 2 (R) ; the operator f(H) 
belongs to C^(A). 
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2.4 The trace per unit volume on A 

Given a T-invariant probability measure P on Q, a positive trace 7b on A (and each A B ) is defined 



The following is readily verified. 

Lemma 1 7b is a linear functional defined on all of A and satisfies 

(i) (normalization) 7b (1) = L 

(ii) (positivity) 7b (^4*^4) > and T B (A*) = T B (A) 
(hi) (cyclicity) T B (AB) = T B (BA) 

(iv) (norm bound) T B (\AB\) < \\A\\T B (\B\) where \B\ = (B*B)^ 

(v) (invariance) T B (V A) = for A G C*(A B ) 

(vi) (partial integration) T B (AVB) = -T B (S7AB) for A,B G C^(A B ) 

If P is in addition ergodic, then Birkhoff 's ergodic theorem implies that for any increasing sequence 
(A/)/ e N of cubes centered at the origin 



for P-almost all uj G Vt. This shows that 7b is the trace per unit volume. 
2.5 The derivative w.r.t. the magnetic field 

For sake of concreteness, we restrict to the cases d = 2 and d = 3 so that B is parametrized by Bj, 
j = 1, 2, 3, as given in ([8]). It is straightforward to treat higher dimensions, but as it leads to heavier 
notations we refrain from doing so here. On C^(H x Q x Z d , Mat(L, C)), the set of all continuously 
differentiable matrix-valued functions with compact support, let us set for each j = 1, 2, 3, 



We also write 6 = (Si, ... ,60)- The operator 6j is called Ito derivative in |Bel2t IBel3] . It is not a 
derivation, but rather satisfies the following algebraic identities. 

Lemma 2 Let d = 2 or d = 3. Whenever all objects below are well-defined, one has: 
(i) With the index j taken mod 3, 



by 




(18) 




(S,A)(B) = d Bj A(B). 
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and 

&j [A, B] = [6jA,B] + [A, SjB) + ~(\y j+1 A,V j+2 B] - [V j+2 A,V j+1 B]) . 

(ii) For an invertible A £ A, 

5,{A- 1 ) = A' 1 [~S S A + % - V^AA^V^A - % - V ) +2 A A^V ' j+1 Aj A' 1 . 

{si) W) = 6MY 

(iv) 5jV k A = V k 5jA 

(v) [generalized DuHamel formula) For z G C and index j taken mod 3, 

5 je zA = z I dse^-'^djAe 9 ** 
Jo 

+ l -z 2 [ ds [ dr e ^~ s > A {V j+1 A e rzA V ' j+2 A e ( s ~ r > A - V j+2 A e rzA V j+1 A e {s ~ r)zA ) 
2 Jo Jo 

(vi) d B .r B {A) = r B (^(A)) 

Proof, (i) Let us focus on j = 3. Deriving ( I14jl using fl9]) shows 
(6 3 {AB) - (S 3 A)B - A(5 3 B)) (B, u,n) = ^ A(B, u, I) B(B, T~ l u, n-l) % - {n x l 2 - n 2 l x ) e^ m . 

Thus replacing n x l 2 — n 2 l\ = (ni — l\)l 2 — (n 2 — l 2 )l\ and then the definition of Vi and V2 allows 
to complete the calculation leading to the first identity in (i), which implies the second one. Then 
(ii) follows from 5j{A~ 1 A) = 0. The points (iii), (iv) and also (vi) follow immediately from the 
definitions. The proof of (v) mimics that of DuHamel's formula, e.g. [Sak[ p. 68]. For n > 1, 

5,e zA = 5 3 {e- A ) n 

= 6 i ei A (eZ A ) nr - 1 + eS A 6 j {ei A ) n - 1 + Jv i+ ie^V j+2 (e5 A ) B " 1 - \ V j+2 e^ A V i+1 (e* A ) B - 1 

n-1 

= e^ zA (Sj et A e^ A + \ V, +1 ei A V J+2 e^ A - \ V, +2 V i+1 e^ A ) , (20) 

k=0 

where the third identity follows from iterating the second one. Now, with norm convergence, 

lim n5je~ A = z5jA, lim nVje™" 4 = zVjA . 

n— >oo n— ¥00 

Therefore the norm convergent Riemann integrals in (1201) lead to 

5j e zA = [ ds e (1 - s)zA (z 5jA e szA + -z V J+1 A V j+2 e szA --z V j+2 A V j+l e szA ) . 
Jo V 2 2 / 
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Finally replacing the remaining derivations of exponentials by their DuHamel formula and changing 
variables readily allows to complete the proof. □ 

As an illustration, let us point out that for H = (H^^n in the form (ITUl) . one has 5jH = 
because ( 1T51) is independent of B. As the products ( 1T4|) explicitly depend on B, this does not imply 
SjH 2 = though. In fact, SjH 2 = ^\Vj+iH, V3+2H]. Derivatives of higher powers can be calculated 
iteratively and expressed in terms of Vj+iH, and functions of H. 

Next let us furnish C£(S x x Z d , Mat(L, C)) with the norm 

1Mb = \\A\\ +^||V^|| 

3=1 3=1 

which dominates the norm in C^(A). The closure w.r.t. this norm is denoted by C 1 ( v 4.). Lemma [2]^i) 
allows to check ||AB||cn < ||A||cn||.B||ci. As also ||^4.* He? 1 = H^llc 1 ? C l (A) is actually a Banach *- 
algebra. In particular, C X (A) is invariant under holomorphic functional calculus. Thinking of Banach 
algebras of analytic functions, invariance under holomorphic functional calculus is the optimal result 
of general nature that one can expect. Here one, moreover, has the following fact. 

Proposition 2 Let H = H* be an element of C 1 (A). Then for any compactly supported smooth 
function f , the operator f(H) belongs to C 1 (A). 

Proof. Let us provide two arguments. The first one is based on the Helffer-Sjorstrand formula as 
presented by Davies [Davj . In fact, Lemma [2]^ii) implies 

WSj (z - H)-^ < C\%m{z)\- 3 , 

so that hypothesis (H3) of [Davj holds, which then directly implies the result. For the second proof, 
let us use that Lemma [2](ii) shows 

\\6 je ltH \\ < C(t + t 2 ) . (21) 

Moreover, one has || Vj e ttH \\ < Ct so that the bound fT2~TT) actually also holds in the C 1 -norm. Now 
one has 

f(H) = Jdt f{t) e* H , fit) = e~* Et fiE) , (22) 

where the integral is norm convergent in A because / is a Schwartz function. Therefore the above 
estimates show that fiH) is in C 1 ( v 4.). □ 

2.6 The von Neumann algebra and the Sobolev spaces 

For p G [l,oo), the Banach space £ p (.Ab,7b) is the closure of Ab under the norm ||^4||lp = 
(r B (|A|f)) 1/p . if 7r GNS denotes the GNS representation of T on L 2 (^4b, 7b), -^ 00 (^b ) '7b) denotes 
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the von Neumann algebra 7t gns (^4b) // where " is the bicommutant. By a theorem of Connes [Con], 
L°°(Ab, 7b) is canonically isomorphic to the von Neumann algebra of P-essentially bounded, weakly 
measurable and covariant families (A^^n of operators on "H endowed with the norm 

\\A\\ L oo = P-esssup HAJbch) . 

Consequently, the family of representations ir^ extends as a family of weakly measurable representa- 
tions of L°°(^4b? 7b)- Moreover, the trace 7b extends to L°°(^4b? 7b) as a normalized normal trace. 
Then L 1 (^4b? 7b) can be identified with the predual of L°°(^4b, 7b), namely the linear space spanned 
by the set of normal states on L°°(^4b, 7b)- Furthermore, one disposes of an associated Sobolev 
space if 1 (,4b? 7b) which is given by the closure of As w.r.t. the norm 

II^Hfli = (r B (|A| 2 ) + r B (|vA| 2 ))^ . (23) 

It is a two-sided ideal in L°°(Ab, 7b)- The Sobolev space is closely linked to dynamical localization, 
see Section [2T91 below and [BESJ. 

2.7 Time evolution 

Let us first consider the Hamiltonian time evolution associated to a fixed self-adjoint element H G Ab- 
Then the associated Liouville operator Ch acting on Ab is defined by 

C H (A) = i[H,A] . (24) 

It is a bounded *-derivation of Ab generating a one-parameter group of ^-automorphisms telK 
e H in Ab- Then observables satisfy the Heisenberg equation dtA(t) = i[H, A(t)]. There is a dual 
time-evolution of states on Ab- Let us focus on normal states w.r.t. 7b represented by positive 
elements p G Ab satisfying the normalization condition 7b (p) = 1- The state is then given by 
A G .4b l— 7b(p^4)- Now the Hamiltonian time evolution of such states is governed by the Liouville 
equation d t p(t) = i\p(t),H], namely p(t) = e~ tCl1 (p(0)). 

For the theory of orbital polarization, one needs also to consider time-evolutions governed by a 
time-dependent Hamiltonian t G [0,T] t— > H(t) given by a continuous path of self-adjoint elements 
in Ab- Then obeservables in Ab evolve according to the time-dependent Heisenberg equation 

d t A(t) = C m (A{t)), C H{t) {A) = i[H(t),A], 

and normal states according to the Liouville equation 

d t p{t) = -£ m (p(t)). (25) 

As to existence and uniqueness of these equations, let us focus on f l25|) . Because H(t) is self- 
adjoint, the Liouville operator C-n(t) is a bounded ^-derivation. Using the Dyson series |RS2 
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Theorem X.69], one readily shows that C-H(t) generates an evolution family (t, s) G [0, T} 2 h-> r] t s of 
*-automorphisms satisfying r\ t t = 1 and r] tjr r] rs = r] t s for s,r,t G [0, T]. The unique solution of (I2l)j) in 
C 1 ([0,T],.A) with p(0) = p ^ A is given by = ?7t,o(Po)- As the generators are bounded, the map 
(t, s) G [0,T] 2 i — ^ r] t s (A) G is differentiable. If one assumes in addition that H(t) G C]{A) for all 
t G [0, T], then are bounded operators in the Banach space C^(A). Hence, if the map t \-t H(t) 
is also continuous w.r.t. the norm in C^(A), then, again by inspecting the Dyson series, one sees that 
the evolution family r] ttS leaves C}{A) invariant and that also (t, s) G [0,T] 2 i-> T)t }S (A) G C^(A) is 
differentiable. Two facts about the solutions of ( 1251) are now of physical importance. First of all, as r) t , s 
is a *-automorphism it follows that the particle density T(p(t)) = T(^ i o(po)) — 7~(Po) is conserved. 
Second of all, if p = Po is a projection, then also p{t) 2 = r]t,o(po)Vt,o(po) = Vt,o(Po) = Vt,o{Po) = 

2.8 Current operator and current-current correlation measure 

Pysically, the current is the derivation of the position w.r.t. time. As these are all observables 
and time derivatives are thus given by Heisenberg's equation, it follows that X = i[ti-b^{H), X] = 
7TB,o;(Vif) where the second equality follows from ffTTj) . Hence the current operator J is defined for 
any H G ^(^4) as an element of A: 

J = WH . (26) 

The following result implies, in particular, that the current vanishes in all states which are functions 
of H. 

Proposition 3 Let H = H* be in C^(A) and f,g be the restrictions to the spectrum of H of 
continuously differentiable function on R. Then 

T B (g{H)Vf(H)) = . (27) 

Proof. As the trace 7b is invariant under V, it follows from Leibniz rule that = 7b(V-£P) = 
nTB(H n ~ l V 'H) for all n > 1. Thus by density 7B(/(-f^)Vi7) = for any continuous function. From 
this follows in turn that 7b(/ {H)V H m ) = for all m > 0. Again by density the result follows. □ 

Next let us recall |KPl [SBBJ the definition and a few properties of the current-current correlation 
measure m. Whenever the current operators VH are in the algebra, or otherwise stated H G C^(A), 
then for each B this measure is given by 

T B (f(H)VHg(H)-VH) = [ m(dE,dE') f(E) g(E') , f,geC (R), 

JRxR 

where the dot denotes the scalar product VHg(H) ■ VH = Yuj=i ^ 3^g{H) VjH. Therefore m is a 
positive measure on M 2 with finite mass Tb{WH ■ VH). Now by Cauchy-Schwarz 

|7b (f(H)ViHg(H) V 3 H) | 2 < T B (|if(#) \g\ 2 (H) V»#) T B (V,H V ,H) 
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and thus, by the Radon-Nykodym theorem, there exist functions e L 1 (R 2 ,m) such that 
T B (f(H)ViHg{H) VjH) = f m(dE, dE') D id (E, E') f(E) g(E') . 

JRxR 

Let us set m^j = D^m. Only the diagonal measures rrijj are positive, but the matrix- valued measure 
( m i,j)i,j=i,...,d is again positive. Furthermore, if A x A' is a rectangle, then again by Cauchy-Schwarz 

m(AxA') < (Tn({P A VH)*P A VH)Tv({P A 'VH)*P AI VH)y < (7b(Pa) MP A >)f 2 \\VH\\ 2 . 
Therefore, if the density of states measure M is defined by 

T B (/(#)) = f M{dE) f(E) , / G C (R) , 
Jr 

then 

m(AxA') < C (JV(A)JV(A'))3 . 

This shows ct-Holder continuity of m if jV is 2a-H61der continuous. The current-current measure 
allows to write out formulas for the conductivity in the relaxation time approximation and its behavior 
on the diagonal E — E' is linked to dynamical properties of H which allows to prove the anomalous 
Drude formula [SBB] , Just below in Section I2l3| the link of m with dynamical localization will be 
explained. 



2.9 Dynamical localization 

This section briefly resumes the notion of dynamical Anderson localization. This is a property of 
random Hamiltonians ttb,uj{H) at a fixed magnetic field B. While there are numerous characteriza- 
tions of this localized phase, here we will merely use the notion of dynamical localization in spectral 
(Borel) subset A C R in the sense that the localization length defined by 

f T dt 

f(A) = sup / -ETr (0| X A(7rB,.(^))(e^ Wi Xe-^^-X) 2 XA (7TB,.(^))|0) 

T>0 Jo 1 

is finite. For quantum Hall systems, this is a sufficient condition for the existence of plateaux [BESJ 
and it will also emerge below as natural condition for the study of magnetization. No attempt is made 
here to prove it in concrete models, but it is known to hold in a strong disorder or band edge regime 
(see jBES] for some concrete verification using the fractional moment method). It will be useful 
that the dynamical localization condition £ 2 (A) < oo can also be checked on the current-current 
correlation measure, namely one has 

£ 2 (A) = 2 f m(dE,dE') 1 
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Now by [BESt Section 5.2], £ 2 (A) < oo implies that the Fermi projection P M = /oo )At (i7) is in 
H l (A&, 7b) for every /x G A (and depends continuously on //). Furthermore, there exists a A/"- 
integrable function E G A i— > £ 2 (E) such that 



2.10 The time reversal operator 

Time reversal is implemented in A by the involutive anti-automorphism 77 : A — > A defined by 



where I = e msV is the rotation is spin space by 180 degree (which is unitary) and the overline denotes 
the complex conjugation. Note that rf = 1, rj(XA) = \r](A) and I 2 = ±1 pending on whether the 
spin is even or odd. Thus one has 



Mv(A)) = T-b(A) , V (VA) = -V(v(A)) , V (SA) = -S( V (A)) . (28) 

3 Orbital polarization 



Following up on the discussion in the introduction, this section presents the main mathematical results 
on the orbital polarization. Let be given a time-dependent Hamiltonian t G [0,T] h-> H(t) G C^(A) 
at fixed magnetic field B, namely a continuous path of self-adjoint elements in Ab- Then states 
evolve according to the Liouville equation (1251) . The change AP in polarization during a change of 
the system between time t = and t — T is now defined by (pQ). It follows from Propoisition [3] that 
the average current density 7b (p J) vanishes for any equilibrium state p = f(H). However, changes 
of the Hamiltonian in time can give rise to non-vanishing currents so that AP does not necessarily 
vanish. Our first result is a formula for the current density solely in terms of the time-evolved state 
p(t), assuming that p(0) = P is a projection in C^(A). 

Proposition 4 Let t i-> H(t) be a continuous path in C^(A) and let P G C^(A) be any projection. 
Then the current associated to the solution P(t) = rjtfl(P) of the Liouville equation ( |25i) can be 
written as 



Proof. Recall that P(t) G C*(A) for all t, so the expression on the right hand side makes sense. For 
sake of notational simplicity, we suppress all t-dependencies in the following computation. Starting 
from the right hand side, one finds 




for all Borel sets A'c A. 



( V {A))(B,cj,n) 



I A(— B, u,n) I , 



T B (P(t)VH(t)) = iT B (P(t)[d t P(t),VP(t)}) . 



*T B (P [d t P, VP]) 



T B (P[P,P]VP) - 7b (P (VP) [H, P]) 
- 7b ((VP) PH) - 7b ((VP) HP) 
T B {PV{PH)) - T b (P(VP)P0 
Tb(PVP) . 
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Here the first equality uses that P(t) solves the Liouville equation (|25|) , the second equality follows 
from the identity P(VP)P = that holds for any projection P, the third one is integration by parts 
from Lemma [T] (vi) and the final one is the Leibniz rule. □ 

Physically the relevant initial state is the Fermi projection. By Proposition [T] all spectral projec- 
tions of H that belong to a part of the spectrum separated by a gap from the rest belong to C^(A). 
Hence, Proposition H] is applicable if initially the Fermi level p lies in a gap of the spectrum of H(0). 

Corollary 1 Assume that the Fermi level p lies in a gap of the spectrum of H(0) and let po = 
x{H{0) — p), where \ denotes the characteristic function on the negative numbers. Then 

AP = i [ dtT B (p(t)[d t p(t),Vp(t)]) . (29) 
Jo 

The formula (|2"91 holds quite generally, but it requires knowledge of the solution of the Liouville 
equation. In order to express the polarization AP in terms of the spectral projections of H(t), let 
us consider (1251) in the adiabatic limit, namely the Liouville equation 

ed t p(t) = i\p(t),H(t)] , (30) 

with a small adiabatic parameter e < 1. If one further supposes that the gap at the initial Fermi 
level /i(0) remains open for all times t G [0,T], i.e. that there is a continuous function p(t) with 
p(t) (fi a(H(t)) for all t G [0,T], then the adiabatic theorem implies 

pit) = P (t) + 0(e) , P (t) = X (H(t)-p(t)) . 

At least formally, this immediately yields a generalization of the King-Smith and Vanderbilt for- 
mula ([3]) with errors of order 0(e). However, using superadiabatic approximations and a homotopy 
argument, a much stronger statement can actually be proved. Note that plugging the adiabatic 
approximation directly into the formula for the current yields the wrong result, since in the presence 
of a gap Po(t) = f(H(t)) for a suitable smooth continuation / of the characteristic function and thus, 
according to Proposition [3], one has T(Poit)'VH) = 0. 

Theorem 1 Let 

(i) either H e C N+2 ([0,T},C^(A)) with 

d? H(t)\ t =o = d? H{t)\ t=T = , for all 1 < n < N , 

(ii) orHe C N+2 (R,C^(A)) be T -periodic. 

Further suppose that there is a continuous function p(t) with p(t) ^ a(H(t)) for all t G [0,T]. Let 
Po = x(H(0)— /i(0)) in case (i) and p = Pn(Q) i n case (h), where P^(0) is a superadiabatic projection 
constructed in the proof. Then 

AP = i f dt T B (Po(t) [d t P {t),VP {t)\) + 0(e N ) . (31) 
Jo 
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Proof. Theorem [9] in the appendix provides superadiabatic projections P^(t) which according to 
Corollary approximate the solution p(t) of (1301) up to small errors in the sense that 

\\P £ N (t) ~ P(t)\\ + ||V(J$(f)-p(f))|| + \\d t (P^(t) - p(t))\\ = 0(e N ). 
With and Lemma IU(iv) this yields that 

AP = z f dtT B {P £ N {t)[d t P £ N {t),VP £ N {t)}) +0{e N ) . 
Jo 

But by the following calculation, the above integral is independent of e so that one can pass to the 
limit e — > using (|53|) without changing its value. For case (i) let us use that the vanishing of the 
first N derivatives of (4 H(t) at the endpoints implies P^(0) = P (0) and P£(T) = P (T). With 
the notation P^ = d t P^, one finds 

<9 £ £ dt T B (p £ n [P £ n , VP£]) = £ dt T B (p£[9 6 p£, VP £ n ] + P £ N [P £ N , d £ VP £ N }) 

= r B (p £ N [d £ p £ N ,vp e N ])\ T Q - J dtr B (p £ N [d £ p £ N ,vP e N ])+ J dtr B (p £ N [p £ N ,vd £ p £ N ]) = o. 

In both equalities it was used that all the derivatives P^, d e Pfj and VP^r of the projection Pfj are 
off-diagonal with respect to P^ and thus l B (d £ P^P^ [ 'V P^) = 0. In the second equality we used that 
differentiability of the map e >->■ P^ from [0, £jv) to C^QO, T], Cg(A)) implies existence and equality 
of the mixed derivatives, d £ P^ = d t d £ P^ and <9 e VP^ = Vc\P^. The last equality follows because 
the boundary terms vanish or cancel each other (in cases (i) and (ii) respectively) and the integrands 
in the last two integrals cancel each other as shows a partial integration w.r.t. V. □ 

Furthermore, the leading order term of (l3~Tj) is invariant under diffeotopies. A differentiate path 
Pi : [0, T] — > Cg(A) of projections is called diffeotopic to Po whenever there exists a differentiable 
map F : [0, 1] -> C7 X ([0, T], C*{A)) such that F(0,t) = P (t) and F(l,t) = P^t). Replacing e in the 
proof of Theorem [T] by the diffeotopy parameter, one obtains the following. 

Corollary 2 Let either of the hypothesis in Theorem [1] hold and suppose that the diffeotopy between 
the paths P and P\ satisfies, in case (i) F(a, 0) = Po(0) and F(a,T) = Po(T) for all a £ [0, 1], or, 
m case (ii), P(a,0) = F(a,T) for alia G [0,1]. Then 

dtr B (Po{t)[d t P {t),VPo{t)})= [ a , tr B (Pi(t)[9iPi(t),VPi(t)]). 

This explains the experimental fact that the total polarization depends only on the initial and 
final configuration of the solid. However, this only holds if two paths are diffeotopic. For general 
cyclic deformations, i.e. H(T) = H(0) and thus Pq{T) = Po(0), the super-adiabatic approximation 
of the charge transported in one cycle 

[ T dtr B (P (t)[d t P (t),VP (t)]) G 2nZ d 
Jo 
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is actually a topological quantity, as was already emphasized by Thouless |Tho] in a more restricted 
context. To see this, let us consider the C*-algebra A = C(Sy) <8> A of continuous functions on 
the circle = [0, T) with values in A. Operators therein will be denoted by A = (A t ) te gi with 
A(t) G A. Besides the derivations V, this algebra has the densely defined ^-derivation id t and a 
trace 

%{A) = [ dtr B {A{t)) . 

Jo 

Thus the polarization for a T-periodic Hamiltonian is given by 

AP = 2 7rCh(P ) + 0(e N ) , 

where P = (Po(^))te§i, e -A an d the vector of associated Chern numbers of a differentiable P G A 
defined by 

Ch(P) = -^f B (P[z^P,VP]) 

Z7T2 

has integer components. This shows that even in the presence of disorder the charge transported 
under cyclic deformations is quantized. 

Corollary 3 Let H\(t) = H (t) + XV(t) be such that H\(t) satisfies the assumptions of Theorem [TJ 
case (ii), for A G [0, 1] with a continuous choice for /j,(t, A). Then the charge transported in a periodic 
cycle of H\(t) is quantized and agrees with that of Ho(t), up to errors of order e N . 

Proof. It follows from standard perturbation arguments, c.f. Lemma [3] in the appendix, that the 
Fermi projections of H\(t) and Ho(t) are diffeotopic in the sense of Theorem [TJ □ 

4 Magnetization at finite temperature 

In this section, the magnetization at finite temperature will be calculated for Hamiltonians in dimen- 
sion d = 2 or d = 3. It will be expressed in terms of the current-current correlation measure. Let us 
begin by recalling the definitions (jl]) and ([5]) of the pressure p((3,fi) and magnetization M(/3,/i) at 
inverse temperature (3 and chemical potential fi. If B = and the time reversal symmetry relation 
f](H) = H holds at B = 0, then one deduces from Lemma Wiy) arid the identities (|28|) that 

M(P,n) = To(<51n(l + e^W-^)) = 7^(ln(l + e^^))) = — M((3, fx) . 

As the magnetization is real, it follows that it vanishes for any time reversal symmetric system so 
that we suppose from now on time reversal symmetry to be broken. This does not necessarily mean 
that B 7^ though as there may be local magnetic fields breaking the time reversal symmetry. In 
order to calculate the magnetization, let us use the identity 

ln(l + e~ pE ) = ln(2) + / d/3' dp, ln(l + e~ p ' E ) = ln(2) + / dp' (1 + e^)" 1 (-E) . 

Jo Jo 
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Thus, if fp tfl (H) = (1 + e^ u ^) 1 is the Fermi-Dirac distribution, the cyclicity of the trace shows 

M,(/3, fi) = i J dp d Bj T B (fp>AH) (fi - H)) . (32) 
Thus it is convenient to set 

Mj(/3,fi) = d Bs T*(ffo(H)Qi-H)) . 

The next result follows from the observation that Mj((3, fi) can be calculated in terms of the current- 
current correlation measure m defined in Section 12.81 

Theorem 2 Let 5H = 0. Then, with j taken cyclically, 

Mj(f3,fi) = % [ m j+w (dE,dE')g^(E,E') , (33) 
Jtsl 2 

where gp jfl is the following smooth function 

1 1 / 1 hid + e -^ E '~^) —hi(l + e-Pl E -ri)\ 
mJE,E>) = + ^^±£ E'-E "J - iE ^ Er (34) 

Proof. In order to simplify the formulas, let us shift H to H — fi so that we may assume that 
/i = and suppress the /i-dependence in the following. Let us also choose a smooth compactly 
supported function / that is equal to fp on the spectrum of H. Its Fourier transform with choice of 
normalization as in (|22|) is denoted by /. Let us also set F = f(H). Then one has 



Mj(f3) = -Tb(H6jF) = - J dt f(t)T B (H5 je ' 



Ht\ 



Now replacing the generalized DuHamel formula, using the cyclicity of the trace and replacing the 
current-current correlation measure shows 

Mj(p) = - j dt f{t) J ds dr % - (it) 2 Tb (H [e* m ^V j+l H, e lHtr V J+2 H}) 

= J dtf{t) J\r(l-r) l -t 2 Jm j+1 , j+2 (dE,dE') ( Ee lEt ^e lE ' tr - E ' e lEH ^ e tEtr ) 
= J dt f{t) l -t 2 J m j+lj+2 (dE, dE') 



\ J m j+1J+2 (dE,dE') 



Ee lEt + E'e %m [E + E') (e lEt - e lE>t ) 
i{E - E')t + (E - E')H 2 
EfjE) + E'f(E') _ (E + E')(f(E')-f{E)) 
E'-E (E' - Ef 



Now let us replace Mj((3) into ( )33|) and evaluate the integral over p. For that purpose one can again 
use f = fp- The following identities are then useful: 

J P dp fp:{E) = - I ln(l + e-f> E ) , j* dp f'^E) = ± ln(l + e^ E ) + | fp{E) . 
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After a short algebraic calculation, one finds 



M(R\ l ! MP Irt ( f(E) + f(E>) 2 ln(l + e-^)-ln(l + e-^) \ 

^(Z 3 ) = 2 _/ m j+1 , j+2 (dE,dE) ^ ^— + — - _ — I • (35) 

Regrouping terms and shifting energy back by fi allows to conclude the proof. The smoothness of 
gp,iii in particular, on the diagonal E — E' is readily checked. □ 

As an illustration, let us consider the case of a 1-periodic operator at zero magnetic field. Then 
Bloch theory applies and „4 B =o = C(T d ) <g> Mat(L x L, C). Moreover, Tb=o(A) = J Jd ^% Tr(A(k)) 
and VjA(k) = djA(k) for a differentiable A where dj denotes the partial derivative w.r.t. kj. The 
Hamiltonian H e C^T*) ® Mat(L x L, C) is given by the matrix- valued function k e T d i— )■ H(k). 
Let E\{k) < . . . < Ei(k) be the corresponding eigenvalues, namely the band functions, and let -P;(A;) 
be the eigenprojection corresponding to Ei(k). In this situation, the magnetization can readily be 
evaluated and gives the following formula stated in [XYFNl ISVXNj . The limit — > oo can be taken 
and this then also shows formula (55) in [ResJ. 

Corollary 4 Let H be 1-periodic, 5H = and B = 0. Suppose that all bands are isolated, i.e. 
Ei{k) < E>z(k) < . . . < E^k) for all k E T d . Further let the Rammal- Wilkinson tensor R^ l \k) and 
curvature of the Berry connection VL^{k) for the Ith Bloch band be given by 

R^{k) = l - {Tr(Pi{k) d t P{k) (H(k) - E l {k)) djPtik)) - Tr(P^) d^k) (H(k) - E^k)) d^k)) 

and 



Then the magnetization is given by 



ftg(fc) = i Tr{P l (k)[d l P l (k),d J P l (k)}) . 



(k) . (36) 



Proof. Then the current-current correlation measure is given by 
Jm id (dE,dE')f(E)g(E') = £ J J^L /(£,(*)) g(Eh(k)) d^djE^k) 

Z,n=l,Z^n ^ ' 

This is replaced into the formula ( 1331) for Mj((3,fi). As 

Jim ^(^O = - = o, 
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the diagonal term vanishes and one finds 
% Jl- r Ah 

Mtf,!*) = g E / 7^^W*).^(*))^M^^ ( 37 ) 

l,n=l,l^n^ ^ 71 ' 

This can be further simplified. To make formulas more simple, we suppress the dependence of all 
objects on k in the following formulas. Deriving P\HP n = and P\P n = 0, one obtains the identities 

= PldjHP n + d j P l E n P n + PlE { d 3 P n , = djPtPn + P x d 3 P n , 

from which a straightforward computation shows 

Tr (Pi djH P n diH) = {E l - E n f Tr (P l P n d^) , 
Therefore, in the splitting of gp jlx as in (135|) . the first term in (157)1 simplifies to 

l,n=l,ly^n ' 1=1 

and the second to 

L L 



]T (ln(l + e-^)-ln(l + e-^)) Tr (pd i+1 P,P n 3 i+2 fl) = ^ In(l + e"^«) 



i+2 • 



o / , , V — 1 - J — 1 - // — V, ' 1 '« "J-r^ '/ a / j 

P l,n=l,tyn 1=1 

The sum of these two thus leads to formula (1361) . □ 

The derivation of (136|) in |XYFN] is based on the semi-classical model for periodic solids including 
first order corrections. A rigorous justification of this approach is the content of |ST] . 



5 Magnetization at zero temperature 

In this section we will focus on the case of zero temperature, namely (3 — > oo. Due to (|32|) . one has 

Mj(oo,fi) = lim Mj(f3,fi) , (38) 

0— >oo 

whenever the limit exists^ Even though a stronger result is proved in Theorem [3] below, let us first 
present a calculation of Mj (oo, fi) = Tb((ijl — H)P^) in the presence of a spectral gap at the Fermi 
level //. This serves as motivation for the following and, moreover, the argument is elucidating. 

Proposition 5 Suppose that the chemical potential \x lies in a gap of the spectrum of H. Let P M = 
foo,fi{H) denote the Fermi projection and suppose that 5H = 0. Then, with j taken cyclically, 

M,(oo,/i) = l -TB((»-H)(l-2P fl )[V j+1 P tl ,V j+2 P fl )) (39) 

= - Z -r B (\fx- H\[V j+l P^V j+2 P^)) . 
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Proof. Because under the gap condition P = P M is smooth, compactly supported function of the 
self-adjoint element H in the Banach *-algebra C 1 (^4), it follows from Proposition |2] that P G C 1 (^4). 
This allows us to use the identities of Lemma [2j Let us start from 

M,(oo,/i) = d Bj T B {^-H)P) = T B {5,(^-H)P)). 

Now the product rule for 5j and SjH = shows 

Mj(oo,P) = T B ((tJ.-H)6 j P) + |T B (V, +1 Gu-#)V i+2 P- V j+2 (ji-H)V j+1 P)) . 

But partial integration and the fact that the derivations Vj all commute shows for any two operators 

TniVjAVfB) = -Tn(AVjV r B) = -T B {AV r VjB) = T B {V r AV jB) . (40) 

Hence 

Mj(oo,f3) = T B ((fi-H)5,P) . 

This is the trace of 5jP multiplied by a function of H. Such a quantity can be calculated due to the 
following fact. If one multiplies £jP = 5iP 2 calculated with the product rule of Lemma [2]^i) from the 
left and right by P shows 

PkPP = [VmP ' Vt+2jP] p ■ (41) 

Similarly, 

(1 - P) 5 t P (1 - P) = ~ (1 - P) [V m P, V J+2 P] (1 - P) . (42) 

This is to be compared with the fact that for any derivation V and projection P one has PVPP = 0, 
so the r.h.s. of the last two equations are correction terms to this. Now one can replace these identities 
into 

Mj{oo,/3) = T B ((fi-H)P5 3 PP) + T B ({fJ l -H){l-P)8 j P{l-P)) , 

in order to conclude the proof of (|39|) . from which the second identity follows immediately. □ 

A natural question is now whether and under which hypothesis Mj(j3, //) converges to the expres- 
sion (139|) in the limit /3 — > oo. If such convergence is given, then (!38l) shows that the zero-temperature 
magnitization itself is given by ( J39l) . This will be shown to be true if there is a gap, but also under 
the more general hypothesis that the Fermi level [i lies in a region of dynamically localized states. 
Indeed, the expression ( |39l) continues to be well-defined if P G H 1 (A B ,T B ) so that the Sobolev 
norm fl23|) is finite (by the way, the same holds for Streda's formula and the Hall conductance). But 
finiteness of the Sobolev norm is known to be linked to finiteness of the dynamical localization length 
£ 2 (A) as is recalled in Section [279l If dynamical localization is given, the next result expresses the 
zero-temperature magnetization in terms of the current-current correlation measure. Formally, this 
is simply obtained by taking the limit of the integrand in ( |33l) . 
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Proposition 6 Suppose that 5H = and the density of states has no atom at the Fermi level ji. 
Moreover, let \i lie in an interval A of dynamical localization, namely £ 2 (A) < oo. Then 

lim Mj((3, fi) = - / m j+1 , j+2 (dE,dE) — — . (43) 

/3^oo i j R2 {.b — h,y 

Proof. First of all, let us point out that the integral on the r.h.s. of the claim is well-defined and 
finite due to the localization hypothesis because rriij is dominated by m. Also let us again shift 
energy so that // = 0. Then the integrand gp in formula (133]) for the magnetization has the following 
limit = lim^oo gp: 

n (F p/\ (E + E')( X (E) - x(E>)) 
because de l'Hopital's rule implies 

ln(l + e-^) „ /rnX 

lim >- = -E X (E). 

P^oo jj 

Therefore, one has to show that the integral w.r.t. to nij + ij +2 of the function hp = gp — g^ vanishes 
in the limit f3 — > oo. For this purpose, let us split the integration region R 2 into three disjoint parts 
Ri, R 2 and R3 where 

Rx = {(E, E') G R 2 1 \E\ < /T r '} U {(E, E') G R 2 1 \E'\ < (3~ v } , 
R 2 = {(E, E') G R 2 1 l-E — E'\ < /3~ N } \ R\ , 
R 3 = R 2 \ (Ri U R 2 ) , 

for some < i] < 1 and iV G N. On the two strips making up Ri, the integrand is bounded by 

C 



\h p (E,E')\ < 



(E - E'f 



For (3 sufficiently large, the interval Ip = [—/3 v , (3 v ] is contained in an interval of dynamical local- 
ization. Therefore, using that m is symmetric in the variables (E,E'), 

[ miij {dE,dE') hp{E,E') < 2 I m{dE,dE')\hp{E,E')\ < 2C£ 2 (Ip) = [ J\f{dE) £ 2 (E) . 

Because of the continuity of the density of states measure Af, this vanishes in the limit — > 00 for 
every 77 > 0. In the region R 2 , the limit g^ vanishes identically and hence hp = gp. Using that 
fp{E) = Fp(E) is the ^-derivative of Fp(E) — — -5 ln(l + e~ l3E ), a Taylor expansion of F yields 

, /PPx , SjFgE) + F' p {E + 6)) - 2Fp(E + 5) + 2Fp(E) 

hp(E,E + 5) = — = 0{d). 
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Hence the integrand hp = gp is uniformly of order 0(/3 N ) on R 2 . For the estimate on R3, let us 
rewrite the estimate as 

1 / -p rpi\ ME) - xjE) + fp(E') - x(E') Fp(E) - E X {E) - Fp(E') + E' X {E>) 
li(3{E,E) = + 2 



E'-E ~ [E-E'Y 

and then use that 

\fp{E)- X {E)\ = 0{e-? E ) , \F P {E) - E X {E)\ = 0{Ee^ E ) , 
to conclude that for E, E' e #3, \E\ < \\H\\ and \E'\ < \\H\\ 

\hp(E,E')\ < C f3 2N e-^ . 

Hence on the regions R 2 and R3 the integrand is smaller than any power of 1//3. □ 

Finally, one can show that (1391) is indeed the correct expression for the zero temperature mag- 
netization even in a region of dynamical localization. This confirms that the magnetization is of 
geometric origin. 

Theorem 3 Suppose that the chemical potential \i lies in an interval A with finite localization length 
£ 2 (A) and that the density of states has no atom at \i. Let P^ = foo^(H) denote the Fermi projection 
and suppose that 5H = 0. Then, with j taken cyclically, 

Mj{oo,n) = l -M(»- H)(l-2P,)[V j+1 P fM ,V j+2 P,}) (44) 

= - l -r B {\fi-H\[V j+1 P^V j+2 P,]) . (45) 

Proof. Again we set P = P M . Let I denote the r.h.s. of (|44p . It is split into four contributions: 
/ = ^ T B ((/i - ff)(l - P) V i+1 FP V j+2 ?(1 - P) - (j + 1 + 2)) 
- % - T B ((// - H)PV ]+l P (1 - P) V j+2 PP + + 2)) . 
Focussing on the first one, it is useful to set 

J = T B ((/i-#)(l-P)V J+1 PPV J+2 P(l-P)) . 
Next let Q v = x(H — \i — 2r]) for 77 > 0. Then PQ V = and (l + P)^ = Q rr Moreover, by hypothesis 

lim 7b (1 -P-Q v ) = lim / N{dE) = 0. 
Therefore, using Lemma DJiv) one has 

J = lim TB((li-H)Q v V j+l PPV j+2 PQv 
974,0 V ' / 
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Next let us use 

= Vj[H,P] = [V 3 H,P] - iC H {V 3 P) , 
which implies that for any e > 

(C H + e)(PV 3 PQ v ) = iPV 3 HQ v + ePVjPQr,. (46) 

In order to calculate PVjPQ v , one now wants to invert C H + e. This is indeed possible because Ch 
is an anti-self-adjoint super-operator on L 2 (A, 7b)- Moreover, one has: 

Claim: For any bounded operator B, the equation Cn{PAQ ri ) = PBQ r] has a unique bounded 
solution PA Q r , = C H X (PB Q v ). 

Indeed, the equation can be written in the form 

P(H — fi — T])(PAQ V ) - (PAQ v )Q ri (H - fi-rj) = % PB Q v . 
Now P(H — \i — t\) < —T] P and Q V (H — ji — rj) > r]Q v . Therefore 

POO 

S= dte tp{H -^ ] {iPBQAe-^ 1 *-^ , 
Jo 

is well-defined because the integral converges. But one readily checks £h(S) = PBQ V so that the 
S = C~^{PBQ n ). Moreover, Ch has a large kernel, namely all functions of H, but the off-diagonal 
part in the direct sum representation of P and 1 — P does not contain any vector in the kernel. This 
proves the claim. 

Now inverting the operator Ch + e in fj4*6l) and taking the limit e -)• 0, one finds 

PV 3 PQ, n = i(C H y l {PVjH Q r] ) . 

Similarly, 

Q V W 3 PP = -tiCnr'iQ^jHP) . 
Now let us replace these identities into the formula for J: 

J = lim rJ^-H){C H )-\Q,V 3 HP){C H y\PV 3 HQ r ^) . 
Dealing in a similar manner with the three other terms, one thus finds 

I = lim % - T B ((fi - H)Q V C H \V 3+l H)PC H \V 3+2 H) Q v + + 2)) 
774,0 z v / 

- ~ T B ((/i - H)P C H \V 3+l H) Q n C H \V 3+1 H) P - (j + 1 o j + 2)) . 

As there are now gradients of H, one can directly write out / using the current- current correla- 
tion measure. The Liouville operator gives energy differences E — E', its inverse inverses of these 
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differences. Therefore, upon taking also the limit 77 — >■ 0, 

* f , JC , JE ^ (ji- E) X (ji- fyxiE 1 - fi) - (ji- ETjxfa- E'ME - fi) 



I = % - f m j+lij+2 (dE,dE') 



(£' - P) 2 

* /■ , JC , JE ,n (jjl- E) X {E - n) X (jJL- E') - (ji- EWE* - (i) X (ji- E) 



/ m j+ i >j+2 (dE,dE 



2 y R3 "-j™™»->»~> {E'-EY 

Regrouping and evaluating terms shows that J is equal to the r.h.s. of ( )43l) so that Proposition [6] 
concludes the proof. □ 

Remark Let us show that formula (14"4"I) reduces to the formula found in [CTVRj , more precisely the 
equations (44) and (45) in this work with H replaced by H — /i (these authors also argue that this 
is necessary further below in their work). Actually, |CTVR] only considers Bloch electrons and thus 
formulas in the quasi-momentum representation, but if one transposes their formulas (44) and (45) 
combined with (34) and (35) read in our notations 

M,(cx),/i) = -^7i(v i+ iP M Q M (^-^)v i+2 p^-v i+2 p M g M (//-iy)v j+2 pJ 

- % - T B (Gu - H)V J+1 P, Q,V J+2 P, -(ji- H)V J+2 P, Q M V i+1 Pj , 

where = 1 — P^. Using VjP^Q^ = P^VjP^ twice in the second summand and regrouping the 
terms then directly leads to (1351) . o 



6 Derivative of magnetization w.r.t. chemical potential 

It follows immediately from Theorem [2] or formula (|35|) that 

B, M i(M = '- f mi+1 , j+2 (iE,iE-) {M^Jk^l + 2 fa(g-fa.(g) j . (47) 

From this equation the following link of the magnetization with topological invariants can be estab- 
lished. 

Theorem 4 Suppose that the chemical potential \x lies in an interval A with finite localization length 
£ 2 (A) and that the density of states has no atom at fi. Let P^ = foo,n(H) denote the Fermi projection 
and suppose that SH = 0. Then, with j taken cyclically, 

<^M,(oo,/i) = z7B(P M [V i+ iP M ,V i+2 P M ]) • (48) 

Proof. One has to take the limit of (|47|) which has to be done with some care as in the proof 
of Proposition [6j In particular, the first summand in (T4"T]) vanishes in the limit (3 — » 00 (due to 
dynamical localization and Holder continuity of the density of states). Therefore 



<9 M Mj(oo, /i) = 1 J mj + i j +2 (dE,dE 



„ xfa -E)- x{» - E') 



(E' - E) 2 
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Now one repeat the argument in the proof of Theorem [3] in order to complete the proof. □ 

On first sight, equation ( 1481) is puzzeling because, even if the Fermi level /i lies in a gap of the 
spectrum of H, the r.h.s. and thus the derivative of the magnetization w.r.t. to the Fermi level can 
be non- vanishing. The reason is that the magnetization itself is defined as the rate of change of the 
pressure w.r.t. the magnetic field and the pressure in turn does vary with the chemical potential: 

= Mf(*A H ))- ( 49 ) 

Thus also the magnetiszation varies with the chemical potential (also at zero temperature as shown 
in Theorem [7] below). Moreover, these derivatives are quantized as will be discussed next. Indeed, 
the r.h.s. of (|48|) can be expressed in terms of the Chern invariants of the Fermi projection which 
for any P G H l {A^, 7b) are defined by 

Ch,(P) = 2tuTb(P[V j+1 P,V, +2 P]) . (50) 

The following are the fundamental results on these objects. For d = 2 and L — 1, they are proved 
in [BESj . We refrain from giving detailed proofs for the case L > 1 and d = 3 because they are 
completely analogous. 

Theorem 5 For any P G H 1 (A&,Tb) , the Chern invariant Chj(P) is an integer number given by 
the Fredholm index of a certain Fredholm operator. In particular, Chj(P) is a homotopy invariant. 

Theorem 6 Suppose that £ 2 (A) < oo for some interval A. Let P M = foo,^{H) denote the Fermi 
projection. Then fi G A i— >• Ch^P^) is constant. 

If a model has its Fermi level /i in a region of dynamical localization, it is called an insulator. 
Indeed, the direct conductivity vanishes in the zero-dissipation limit jSBBj . If, moreover, the Chern 
invariant Chj(P M ) does not vanish, one speaks of a Chern insulator, which is a subclass of topological 
insulators. The most prominent (prototypical) example is the Haldane model [Hal] . An insulator 
with vanishing Chern invariants is often also called a normal insulator (there are, however, also 
non-trivial topological insulators with vanishing Chern invariants). 

The Hall conductivity Cij((3, //), i ^ j, is the linear response coefficient of the dissipationless trans- 
verse current. At zero temperature and if the Fermi level lies in a region of dynamical localization, 
it is up to numerical constants given by the Chern invariant 

o- jJ+1 (oo,n) = Ch i+2 (P /i ) . 

For d = 2, a proof is contained in [BESj, the three-dimensional generalization will be dealt with 
elsewhere. It therefore follows from pHj) that for fj, in a region of dynamical localization 

dp Mj(oo,fi) = — <T j+1J+2 (oo,n) . 

The Chern invariant is related to another quantity, namely the derivative of the density of states 
w.r.t. the magnetic field. In the context of the quantum Hall effect, this is known as Streda's formula 
[Str ]lBeT2l . 



26 



Theorem 7 Suppose that the Fermi level \i lies in a region of dynamical localization. Let P^ = 
foo,fi(H) denote the Fermi projection and suppose that 5H = 0. Then, with j taken cyclically, 
Streda 's formula holds: 

d Bo T B (P,) = -iT B {P,[V 1+1 P„V J+2 P,]) . (51) 

Proof. Because it is short and elucidating, let us provide a proof only in the case where [i lies 
in a gap of the spectrum of H (the general case goes along the lines above and is again somewhat 
lengthy). As above, with the notation P = P^, 

Ss,MP) = Tb(^) = 2T B (P6 jP ) + iT B ([V 3+1 P, W1 ) = 2 M PS jP ) . 

Then (HI]) already concludes the proof. □ 
As a corollary, for /i in a region of dynamical localization, 

d^Mjioo,^) = -<9 B .T B (P M ) . 

This follows also if one derives the zero-temperature limit of the identity (|49|) w.r.t. the magnetic 
field. Thus, one can derive Streda's formula from Theorem H] and vica versa. 



7 Boundary currents in Chern insulators 

In this brief section, we sketch how non-trivial Chern numbers in Chern insulators lead to edge 
currents, just as they do in quantum Hall systems (actually, quantum Hall systems are special cases 
of Chern insulators). The results are given without detailed proofs and even without a detailed 
description of the set-up. In dimension d = 2, all this is given in jKRSj, while the case d = 3 will be 
dealt with in a future publication. Let us consider the half-spaces Tj = {(n 1; . . . , rid) \ n j — 0}. Then 
Hj tL0 denotes the restriction of H u = ti-b,uj(H) to i 2 (Tj) by imposing Dirichlet boundary conditions. 
The family (-f/j>) w en is still satisfies the covariance relation ( JTBl in the direction i ^ j. These 
operators are the representations of an operator Hj lying in a Toeplitz extension £ of A which is 
constructed in |KRS] . This extension £ in turn contains operators Aj = (Aj <u ) w€ n for which the 
following trace is well-defined and finite: 

nj>0 

where Cj is the jth unit vector so that rijej is the <i-dimensional vector having vanishing components 
except the jth which is equal to nj. Operators Aj for which \Aj\ has finite 7J-trace are called 
7^-traceclass. The following result can be proved using the techniques of |KRSj (details will be 
presented elsewhere). It shows that the boundary currents of any Chern insulator are determined by 
bulk properties, just as for quantum Hall systems. 



27 



Theorem 8 Let H have a gap A and let F be a positive smooth function supported in A which has 
unit integral f dE F(E) = 1. Then F(Hj) is Tj-traceclass and for i ^ j and any ji e A 

^(F(^)V 4 (^)) = Ch i+i (P M ), 
where j + i is calculated modulo d. 

It may be possible to extend this result to Chern insulators with a mobility gap (Fermi level in 
an interval of dynamical localization) if as in [EGSj an adequate limit procedure for the calculation 
of the edge currents is chosen. 

8 Magnetic susceptibility 

There are two alternatives for the calculation of the orbital magnetic susceptibility: either conven- 
tional linear response theory in the magnetic field, or simply derivation of the magnetization w.r.t. 
the magnetic field. In view of the above, it is natural to follow the second approach which has, 
moreover, the advantage of not needing any approximation. The magnetic susceptibility is then the 
tensor defined by 

XjAPil 1 ) = 9 Bi M j{P,v) , i,j = l,...,d . 
Using (|32|) it is natural to first calculate 

XjdP'V) = d Br 9 Bj T B (f/3,^H) (fi-H)) > 

because XjA^i^) * s then obtained as average over the inverse temperature. In principle, there is no 
problem to calculate /-0, because for any smooth and compactly supported function / with 

Fourier transform /: 

Sif{H) = i Jdt /(f) J ds ^ dr e^ tH [V i+1 He* rtH , V i+2 He« s ^ tH ] . (52) 

This leads to some formula for the susceptibility only invoking the gradients VjH, but the algebraic 
complexity is considerable and it is not clear whether that formula is of any help. On the other 
hand, numerous prior works (nicely reviewed in [BCS] ) have shown that any final expression for the 
susceptibility simply contains many terms so that this complexity seems unavoidable. 

A Adiabatic theorem for C*-algebras 

The following is essentially an adaption of the super-adiabatic theorem of Nenciu |Nen] for unitary 
evolutions on Hilbert spaces to the setting of C*-dynamical systems. All algebras and notations 
are as in Sections [2] and [31 but we emphasize that the results and proofs hold for general C*- 
dynamical systems generated by bounded *-derivations. Hence let t e [0,T] \-t H(t) G C^(A) be a 
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time-dependent Hamiltonian (the derivative 5 w.r.t. the magnetic field plays no role here, so B is 
supposed to be fixed). In the remainder of this appendix, it is always assumed that the following 
hypothesis holds. 

Gap condition: Let cr*(t) C <j(H{t)) be a subset of the spectrum of H(t) such that there exist 
continuous functions f± : [0,T] — > K defining intervals I(t) = [/-(£), f+(t)] such that er*(t) C I(t) 
and 

g = inf dist{I(t),a(H(t))\crJt)) 
te[o,T] 

is strictly positive. 



Lemma 3 Assume the gap condition and let P(t) = XtrM)(H(t)) be the spectral projection of H(t) 
associated with a*(t). Then P(t) G C^A) for all t G [0,T]. If H G C N ([0,T),C^(A)), then also 
PeC N ({0,T],C'(A)). 

Proof. Let j(t) C C be a closed curve in the resolvent set of H(t) encircling <7*(t) once in the 
positive sense with dist (jit), a (H(t)) \ cr*(t)) < |. Then 



P(t) = I ^(z-H(t))-^Cl(A), 



since according to Proposition [1] one has (z — Hit))^ 1 G C^(A) for all z ^ R. The continuity of /±(t) 
implies that j(t + h) is homotopic to j(t) in the resolvent set of H(t + h) for small enough and 
hence 

P(t + h) = I —{z-H{t + h)y x = I — (z-Hit + h))- 1 



j(t+h) 2m 7 7 (t) 2? ™ 



for \h\ small enough. Therefore 



3^ p(t) = I ^dnz-mr 1 e c^a) 



hit) 

for n <, since t h-> #(t) and thus also t 4 (z - Hit))' 1 is iV times differentiable. The latter fact 
follows from 

dtiz-Hit))- 1 = (z-HW-iHWiz-Hit))- 1 , 

and induction. □ 

For the formulation of the following theorem, let us introduce for any differentiable path t G 
[0,T] H> A(t) G Cl{A) the abbreviation 

IPWIII = \\A{t)\\ + ||ftA(t)|| + \\VA(t)\\. 
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Theorem 9 (Superadiabatic projections) Let H £ C N+2 ([Q,T\,C*(A)) for some N £ No and 
assume the gap condition. Let P(t) be the spectral projection of H(t) associated with the gapped part 
of the spectrum just as in Lemma [3j Then there are En > 0, Cm < oo and orthogonal projections 
Pffif) £ C s 1 (^) such that the map e \— > P N {-) £ C 2 ([0,T];Cg(A)) is differentiable and the following 
properties hold uniformly in t: 

IK(f) Ill < C N e, (53) 
\\\[ied t -H(t),P £ N (t)]\\\ < C N s N+1 . (54) 
Ifd? H(t)\ t=T = for some r £ [0, T] and all n<N, then P* (r) = P(r). 

The proof of Theorem [9] is based on the following explicit construction. 

Proposition 7 There exist unique functions P n £ C N+2 ~ n ([Q, T], Cg(A)), 1 < n < N, such that 



n=0 



for <m < N with P (t) = P(t) satisfies 



pe pe _ pe , _m+l/0 , /n{ c rn+2\ 



(55) 



with G 



m+l 



1 PnPm+l-n an d 



%edt-H,R 



i - 1 m 



I £ 



m+l 



P 



(56) 



Here and in the following we sometimes drop the t-dependence in the notations to make formulas 
shorter. Ifd?H(t)\ t=T = for some r £ [0,T] and all n<N, then P n (r) = for all 1 < n < N. 

Proof. One constructs Pj by induction in m. Let Pq = P, then Pq = Pq so that 



pe pe 



P E and 



e<9 t -Pf(t),P £ (t) = ieP (t) = 0{e). 



Now assume that Pj are already constructed for j = 0, . . . ,m such that f|55|) and ( 156]) hold. One 
needs to show that there is a unique smooth P m+i (t) such that P m+1 = P m + e m+1 P m+ i satisfies 



pe pe 
r m+l r m+l 



pe I m+2^-y 



C(e m+3 ) 



and 



Using f )55|) one finds that 



z £ <9 t -tf(t),P m+1 (t) 



*£ m+2 P m+1 (t). 



(57) 
(58) 



pe pe 
m+l m+l 



+ £ m+1 (P P m+ l + Pm+lPo) + 0(£ m+i ) 
P m + (G m+1 + P P m+1 + Pm+lPo) + 0(£ m+2 ) 
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In order to obtain (157|) up to order e m+2 , P m +i needs to be chosen such that 

Gm+l + -Po-Pm+1 + Pm+lPo = Pm+1 ■ (59) 

Claim 1: (|55|) implies that [P ,G m+ i] = 0. 
Indeed, with P £ - P = 0(e) it holds that 

} + 0{e m+2 ) 

Since [Po, G m+ J does not depend on e, the claim follows. 

Claim 1 implies that fl59|) and thus also ( 157)) holds if and only if the diagonal blocks of P m +i are 
set equal to 

PoPm+lPo — — PoG-m+lPo an d Pq Pm+lPo = Po~G m+ iP^~ , (60) 

where P^~ = 1 — Pq. Note that (|59|) puts no condition on the off-diagonal blocks. To fix the latter 
let us evaluate the l.h.s. of (!58|) using (156)) : 



ted t -H,P m+1 



- e m+1 (iP m - [P,P m+1 ]) + l£ m+2 P m+1 . 
Comparing with (1581) let us require 

*P m - [P, P m+1 ] = 0. (61) 

The following shows that the diagonal blocks of the left hand side in fl6T|) are automatically zero for 
our choice (16"U|) . 

Claim 2: (ES) and ([56]) imply that 



Po P-m Po — ^Po [H, G m+ 1 ] Po and P ± P m P x = -iP^[ H, G m+ i]P . 
Indeed, from 

e m+1 [H, G m+1 ] = e m+1 [H -ied t , G m+1 ] + 0(e m+2 ) 

[H-ied u {P m ) 2 -P m ] + 0{e m+2 ) 
-ie m+1 (p m P + P P m - P m ) + 0{e m+2 ) 

and the fact that [H, G m+ i] + i(P m P Q + PoP m — P m ) is independent of e, it follows that 

[H, G m+ i] + z(P m P + PoP m — P m ) = . 
Taking the diagonal blocks of this equality proves the claim. 
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The off-diagonal blocks of (|6T|) vanish if and only if 

[iP m -[H } P m+1 ],P ] = 0. (62) 
Let z G C be in the resolvent set of H(t) and put R z (t) — (z — H(t))~ l . Then ( 162]) is equivalent to 

R z [iP m - [H, P m+1 ], P ]R Z = . (63) 
Now using [R z , P ] = one checks 

R z [[H,P m+1 ],P )R z = R z [[{H-z),P m+1 ],P )R z = [[P , P m+1 ], R z ) . 
Hence integrating (I6"3"j) along a curve j(t) encircling a*(t) once in the positive sense gives 
f dz • f dz 

<fi Rz[Pmi Po)Rz — o — [Rz, [Poi Pm+l]) = [Pqi [Pqi ^m+l]] = Pq Pm+lPo + PqPiu+iPq ■ 

J~ ( (t) ^ J~f(t) 

This fixes the off-diagonal terms in P m +i- Summarizing, with the choice 

f dz ■ 

Pm+1 = Pq Gm+lPo ~ PqG„ 1+ iPq + ® — R g [P m , Pq\R z , (64) 

/ 7 (t)27T 

one has (ETJ) and (1581). 

It remains to prove the final statement of Proposition [71 First note that fl"(r) = implies 
i? 2 (r) = and thus Po(r) = and Pi(t) = 0. Again one proceeds by induction. Assume that 
d™H(t)\t= T = for some r e R and all n < m + 1 and that P n ( T ) — for all n < m. Then clearly 
G m+ i(r) = by the definition of G m+ i and 



'7(r) 

To see that P m (r) = 0, note that iteration of 

Pm(r) = dt<t —R z [P m . u P }R z 
Jj(t) ^ 



P m+ i(r) = I ^R z (r)[P m (r),Po(r)}R z (r). 

J Mr) ^ 



t=T 



l ^R z (T)[P m ^(r),P (r)}R z (r) 
/ 7 ( T ) 27r 

/ ^ / ^ 2 Wft(T)[P ra - 2 (r) I Po(r)]i? 2 -(r) ) Po(r)]^(r) 

A(r) 2?F J Mr) ^ 



implies that P m+ i(r) = 0, if <9™ +1 Po(r) = 0. This in turn follows form the assumption that 
d?H(t)\ t=T = for all n < m+ 1. □ 
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Proof of Theorem O It remains to turn P N (t) into a projection. Since by flo3]) there is a constant 
such that \\(Pn) 2 — P N \\ < cn£ N+1 , by the spectral mapping theorem 

a(P N ) c [-c^ 7V+1 ,c^ 7V+1 ]U[l-c^ iV+1 ,l + c J ve Ar+1 ] C 



.1 11 U [2 51 



the latter for e < e^r = (4cat) jv+i . Hence one can define for e < En 

dz 



P £ 

N 



z-\\ 



2m 



[* - Pn)- 1 



where the integral is over the circle \z — 1| = | in the positive sense. Then e h-> P n ( 
C 2 ([0,T};C^(A)) is differentiable, since P N and thus also {z — Pn) 1 Taking norms of 

dz 



P £ — P 



2-1 



2m 



z-P N )-\P N -P)(z-P) 



-i 



and its derivative d t and gradient V implies (I53"|) . Similarly, one can use 



[ied t -H } P i 



z-X\=l 



2m 



iedt-H, (z-P\ 



e\-l 
N) 



2-1 



JV+1 



ied t -H,P s N 



[z ~ Pn)- 1 



z-l\=i 



dz 
2^ 



[z - P^PNiz - Pn)- 1 



□ 



to show (1541) . The last claim follows directly from Proposition [71 
Corollary 5 Let rj e s be the family of *- automorphisms providing the solutions of fl3"01) . Then 

vi (pm) = m) + 

u>i*/i 

|||A £ (t)||| = 0(6*1*1). 

Proof. Let Q £ (t) = ^ (P^(0)) so that A £ (t) = Q £ (t) - P N (t). Then A £ (0) = and A £ (t) satisfies 
the inhomogenous linear equation 

dt A%t) = d t (Q%t)-P N (t)) = ~ lQ £ (t) -P N (t),H(t)} + R(t) = I [A £ (*),tf(t)] + i2(t) , 

with 

#(*) = ^£d t -#(*),P^(*)] . 
Hence by the variation of constants formula 

As |||i2(t)||| < C N e N by Theorem El taking norms then gives |||A £ (t)||| < \t\C N e N . □ 
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